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Recent atmospheric neutrino results from Soudan 2 * 
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for the Soudan 2 Collaboration 

(Argonne National Laboratory, University of Minnesota, Tufts University, USA; Oxford University, 
Rutherford Appleton Laboratory, UK) 

An updated measurement of the atmospheric v^jv^ ratio-of-ratios, 0.68 ± 0.11 ± 0.06, has been obtained using a 
4.6-kty exposure of the Soudan-2 iron tracking calorimeter. The L/E distributions have been analyzed for effects of 
t'p — ► Vx oscillations, and an allowed region in the Am 2 vs. sin 2 2# plane has been determined. 



1. INTRODUCTION 

Soudan 2 is a fine-grained iron tracking calorime- 
ter located 2100 mwe underground in Soudan, Min- 
nesota, USA. The experiment has been taking data 
since 1989 when the detector was one quarter of its 
full size; the construction was completed in 1993. 
The experiment continues to take data with 90% 
live time. The results presented here were obtained 
using a 4.6 fiducial kiloton-year exposure. For pre- 
vious Soudan-2 results on the ratio of ratios of the 
atmospheric neutrinos, 
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seeRefs. [%%. 

Detailed descriptions have been published of the 
design and performance of Soudan-2 calorimeter 
modules [ [| , and of the active shield system [ [| . 

2. ATMOSPHERIC NEUTRINO FLAVOR 
RESULT 

We analyze fully contained events (all hits more 
than 20 cm from the nearest detector edge) to 
determine R using single-track and single-shower 
events for which the active shield registered no in- 
time two-layer coincident hits (gold events). This 
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sample contains mostly quasi-elastic neutrino inter- 
actions, with a background of photon and neutron 
interactions originating in cosmic-ray muon inter- 
actions in the rock. These rock events are mostly 
flagged by hits in the active shield, but some are 
not accompanied by shield hits and constitute a 
background for our neutrino sample. A correction 
for this background is determined using distribu- 
tions of event depth in the detector. For the 4.6- 
kty Soudan 2 exposure, we obtain a ratio-of-ratios 
value of 

R = 0.68 ± 0.11(stat.) ± 0.06{syst.) 

Our full detector simulation of atmospheric neu- 
trino interactions is based on the 1989 Bartol flux 
for the Soudan site [ See Ref. [ |(| for more 
information on the quasi-elastic data analysis. 

3. HIGH RESOLUTION SAMPLE 

We wish to analyze neutrino L/E distributions 
where the neutrino path length is calculated from 
the event zenith angle, radius of the Earth, and 
the neutrino production height. We therefore make 
data selection designed to optimize our resolution 
in event angle and in energy. For the quasi-elastic 
events (tracks, showers) we require the charged lep- 
ton momentum to exceed 150 MeV/c if a recoil pro- 
ton is measured in the event, otherwise we require 
event visible energy to exceed 600 MeV. For multi- 
prongs, we require the visible energy to exceed 700 
MeV, the vector sum of visible momenta to exceed 
450 MeV, and the charged-lepton momentum to 
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exceed 250 MeV. The event counts in the resulting 
event samples are given in Table 1. 

Table 1 

Soudan-2 HiRes event samples for the atmospheric 
neutrino data (before and after the background 
subtraction) and Monte Carlo. 





v„ 


v e 


Data, before 


99.0±10.0 


121.0±11.0 


Data, after 


92.3±10.1 


114.6±11.1 


Monte Carlo 


133.6±4.4 


114.6±4.1 (norm.) 



We calculate our 'resolution' as the difference be- 
tween the Monte Carlo truth and the reconstructed 
final-state characteristics, thereby taking into ac- 
count both the event kinematics and the instru- 
mental resolution. The resulting AE/E is 20%, 
and the average angular resolution is 33° for v^, 
and 21° for v e charged-current events. The corre- 
sponding resolution in \og(L/E) is 0.49 and 0.43, 
respectively. 

In this 'high-resolution' (HiRes) sample, the rock 
background is small, 5-7%, and the neutrino fla- 
vor is determined correctly 92% of the time. The 
Monte Carlo sample represents 25 kty of exposure. 

4. L/E RESULTS 

In Fig. 1 we show the L/E distributions for the 
HiRes sample. Fig. la shows that the ^ e -flavor 
data (crosses) follow the shape of the j/ e -flavor no- 
oscillation Monte Carlo distribution (dashed his- 
togram). (The MC histogram is normalized to the 
v e data.) In the ^-flavor sample, shown in Fig. lb, 
the depletion of the data vis-a-vis the no-oscillation 
Monte Carlo is obvious (with the same MC normal- 
ization as for the v & sample). 

To convert results of our atmospheric neutrino 
simulation generated under the no-oscillation hy- 
pothesis into simulated neutrino oscillation data, 
we apply to every MC event an L/S-dependent 
weight representing the probability of no oscilla- 
tion for a given Am 2 and sin 2 29 valid for two- flavor 
neutrino oscillations. Fig. 2 displays the ^-flavor 
L/E distribution (crosses) together with the oscil- 
lation MC for sin 2 2# = 1 and for four values of 
Am 2 . We see that the simulation exceeds the data 
for Am 2 = 0.0001 eV 2 and 0.001 eV 2 for most of 
the downgoing neutrinos, while the simulation is 
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Figure 1. L/E distributions for Soudan-2 v e - and 
f M -flavor HiRes samples (crosses) together with the 
no-oscillation atmospheric neutrino Monte Carlo 
(dashed histograms). 



smaller than the data for Am 2 = 0.01 eV 2 . The 
best agreement between data and the simulation is 
obtained for Am 2 = 0.008 eV 2 . 

To determine the neutrino oscillation parameters 
Am 2 and sin 2 29 from our data, we fit the atmo- 
spheric neutrino Monte Carlo distribution includ- 
ing the effects of neutrino oscillations to the L/E 
distribution for our data corrected for the rock 
background, by minimizing % 2 . At the same time 
we assume that no neutrino oscillations occur in the 
v e flux, and include only the total v e event counts 
in our % 2 . In addition to Am 2 and sin 2 2#, the 
neutrino flux normalization factor, /„, is the third 
free parameter in the fit. Errors in both the data 
and the Monte Carlo are taken into account when 
calculating x 2 . 

Fig. 3 shows our x 2 surface as a function of 
Am 2 and sin 2 2(9. A distinct minimum is observed 
at Am 2 = 0.008 eV 2 and sin 2 26> = 0.95 (with f u 
— 0.82). The corresponding confidence-level con- 
tours are shown in Fig. 4 for 68% CL and 90% 
CL. The Soudan-2 allowed region is seen to over- 
lap both the Kamiokande and SuperKamiokande 
allowed regions [ . 

Our plans for the near future are to incorpo- 
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Am 2 = 0.0001 eV 2 Am 2 = 0.001 eV 2 




Figure 2. L/E distribution for Soudan-2 ^-flavor 
HiRes sample (crosses) together with the atmo- 
spheric neutrino Monte Carlo with oscillations for 7 - M - Nakahata, These Proceedings. 
sin 2 28 = 1 and for four values of Am 2 (dashed his- 
tograms) . 



rate more recent atmospheric v fluxes in the Monte 
Carlo, and to include partially contained events in 
the analysis. And we certainly plan on taking more 
data, even through the period of dynamite blasting 
of the MINOS cavern in the Soudan mine. We ex- 
pect to reach 5.0 fid. kty by the Millenium. 
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Figure 4. Soudan-2 allowed regions for CL of 68% 
and 90%. 



